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I. Introduction 



Cluster expansions have a reputation of being hard to use; this is largely due to 
the difficulty to capture them in a single short formula. These expansions were introduced 
in constructive theory by Glimm, Jaffe and Spencer [GJS1-2] and they were improved or 
generalized over the years [BF] [BaF] [Bat] [Bl]. For many years the Ecole Polytechnique was 
happy using a cluster "tree formula" due to Brydges, Battle and Federbush (see [Rl-2] and 
references therein). This formula expresses connected amplitudes more naturally as sums 
over "ordered trees" rather than regular trees; but a combinatoric lemma due to Battle and 
Federbush, shows that the sum over all ordered trees corresponding to a given ordinary Cayley 
tree has total weight 1. A slightly disturbing dissymmetry in this formula reveals itself in the 
need to order in an arbitrary way the elements on which the cluster expansion is performed, 
and in the particular role played by the first element or root of the tree. 

Brydges kept digging for a truly satisfying, still more beautiful formula, and with 
Kennedy they obtained it in [BK] (see also [BY][B2]). This formula shows a clear conceptual 
progress; it does not require the use of arbitrary choices such as an arbitrary ordering of the 
objects to decouple, and the outcome can be written in a somewhat shorter form, involving 
directly standard Cayley trees. Both formulas share a positivity preserving property which is 
crucial in constructive theory: the corresponding interpolations of positive matrices remain 
positive. Both can be built by iterating inductively some perturbation step. But the first 
formula insists in completing the cluster containing the "first" cube, before building the next 
one. The second and better formula blindly derives connections, hence all clusters grow 
symmetrically at the same time. Therefore we prefer to call it a forest formula (see below). 

However the original proof of this formula in [BK] relies on a differential equation 
(Hamilton- Jacobi) which is perhaps not totally transparent. The purpose of this paper is to 
derive a fully combinatoric or algebraic proof of this type of formula, and to show how to 
apply it to various examples chosen for their pedagogic value. Many situations in constructive 
theory in fact require several cluster expansions on top of each other, and for this case we 
derive a generalization of the tree or forest formula, which we call the jungle formula. 

We also derive a generalization that applies not only to exponentials of two-body 
interaction potentials, or to perturbations of Gaussian measures. It is a general interpolation 
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formula we call the Taylor forest formula. 

Finally we propose a formula that performs, in a single move, the succession of a 
cluster and a Mayer expansion. 

II. Two forest formulas and their combinatoric proofs 

First let us recall the Brydges-Kennedy formula [BK] under the most convenient 
setting for the following discussion. Let n > 1 , be an integer , I n = {l,...,n}, V n = 
& I n ,i ^ j} (the set of unordered pairs in I n ). Consider n(n — l)/2 elements u^y 
of a commutative Banach algebra £>, indexed by the elements {ij} of V n . An element of V n 
will be called a link, a subset of V n , a graph. A graph # = {Zi, . . . , Z T } containing no loops, 
i.e. no subset {{^i, 12}, {h, 23}, • • • , {ifc, h}} with k > 3 elements, will be called a u-forest 
(unordered forest). A sequence F = . . . , l T ) of links the range of which . . . , l T } is a 
u-forest, will be called an o-forest (ordered forest). 

A u-forest is a union of disconnected trees, the supports of which are disjoint subsets 
of I n called the connected components or clusters of J. Isolated points also form clusters 
reduced to singletons, hence the total number of clusters is n — r, where r = |#|. Two points 
in the same cluster are said connected by 

Now the Brydges-Kennedy forest formula states: 



Theorem II.l [BK] 

\er n ' 3={ii,...,i T } ^u=i Jo ' ^v=i ' ier n 

u — forest 

where the summation extends over all possible lengths r of including r = hence the empty 
forest. To each link of $ is attached a variable of integration hi; and /i^j(h) = inf{/i;,/ e 
L${ij}} where L${ij} is the unique path in the forest $ connecting i to j. If no such path 
exists, by convention h^^h) = 0. 

Our second forest formula writes exactly the same, except that the definition of 
the h^.^h) is different, and it involves for each cluster a particular choice of a root. For 
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simplicity let us slightly restrict this arbitrary choice by imposing the following rule : for each 
non empty subset or cluster C of I n , choose rc, the least element in the natural ordering of 
I n = {1, . . . , n}, to be the root of all the trees with support C that appear in the following 
expansion. Now if i is in some tree X with support C we call the height of i the number of 
links in the unique path of the tree X that goes from i to the root rc- We denote it by l % (i). 
The set of points i with a fixed height k is called the k-th layer of the tree. 
We now have 



Theorem II.2 

ex p( e u ) = e (n f dwi -) (n «0 ex p( e w f ( w )-«o ' ( n - 2 ) 

u — forest 

where the summation extends over all possible lengths r of$, including r = hence the empty 
forest. To each link of$ is attached a variable of integration wi. We define the wf as follows. 
w^ i jy(w) = if i and j are not connected by the If i and j fall in the support C of the 
same tree T of $ then 

w {ij}( w ^ = ® I^W ~~ ^U)\ ^ ^ (i and j in distant layers) 
w^y(w) = 1 if = (i and j in the same layer) 

w {ij}( w ) = w {u'} if - 1 = = l Z (i')> and {H 1 } e X. (i and j in 

neighboring layers, i' is then unique). In particular, if {ij} G then wf..-,(w) = w^jy. 



The proof we give of theorem II. 1 relies on two lemmas. 



Lemma II. 1 Let a,Q, . . . , a p , (p > 1) be distinct complex numbers, then: 



[ 

J t >0,...,t p >0 

t +...+tp=l 



dt . . . dtp exp(t ao + . . . + t p a p ) = ^ 



^ U P j=o{ai -aj) 

i-U ^ 



(II.3) 



Proof: By induction, p = 1 is an easy computation. We assume the result is true for p > 1 
thus 
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/ dt . . . dtp+i exp(t a + • • • + t p+ ia p+1 ) 

J t >0,...,t p + 1 >0 



t +...+t p+1 = 



[ dt p+1 e tp+iap+1 (1 - t p+1 ) p [ dw ... dw p ex.p(w b + ... + w p b p ) (II.4) 

Jq J ivq>0,...,w p >0 



where we performed first the integration on t p+ i and made the following change of variables: 
ti = (1 — t p +i)wi, bi = a,i(l — t p +i) for < % < p. By the induction hypothesis this becomes 



f 1 dt^+^ii-t^r J2 

JO ,_r, 



i=0 



e a * 1 

- u e ap+1 V - (II 5) 

^ Ytftlia, - a 3 ) ^ {UU (a, - a,))(a p+1 - a,) 

This boils down to the wanted expression for p + 1 after remarking that we have the following 
rational fraction decomposition with simple poles: 

i a i i 



X 



(X - a ) . . .{X - a p ) t^(X-di) Yl P i=o (ai ~ aj) 
and putting X = a p +i. 



Lemma II. 2 With the same notation as in the beginning of this section, assume that we 
are given two sets of n(n — l)/2 indeterminates and v^jy; then the following algebraic 
identity is true in the field of rational fractions (D(-U{jj}, v^jy). 

o — forest 

where, forO <v <t, = Y.{ij}/v u {ij}, b„ = {u}/^ meaning that the points 

i and j are connected by the sub-o-forest . . . , l v ) of F . In particular since no points are 
connected by the empty forest a$ = 0, b$ = 0. 



Proof: Both sides of this identity are polynomials in the v/'s, and we prove equality 
coefficient by coefficient. First consider the case of the constant monomial. We must show 

-u h . . .ui T f if n > 2 



£ t- 1 '^ ={ >»= (IL7) 

F=(i 1 ,...,( T ) It (.J 

o — forest 

For n=l it is trivial, so assume n > 2 and denote Af the contribution of F to the left hand 
side. Given an o-forest F = (Z 1? . . . , Z T ), an o-forest F' of the form (Z l7 . . . , / T , / T +i, . . . , l T +n) 
will be called a K-extension of F. When we need not know k we simply say an extension of 
F. 

a F 

Let 7^f = Af • — — where a to tai = J2iev n Ul ^ since n > 2) . 

atotai 

Notice that 

atotai = af + £ m T+1 , (II.8) 

F' 1 — extension of F 

because summing over F' is the same as summing over all links {ij} where i and j lie in 
different clusters of F, while af handles summing over links where i and j are in the same 
cluster. 

Hence 



af = atotai + £ (-^K^ (H.9) 

F' 1-extension of F T + 1 

i.e multiplying by A F /a tota \ 

n F = A F + Yl np " • ( IL1 °) 

F' 1 — extension of F 

Multiple (but finite because r < n — 1) iteration of (11.10) yields 

n F = Af > ( ilu ) 

F' extension of F 

where the sum is over extensions of all possible lengths. In particular 

J2^ F = Yl AF = n$ = A®.— 2- = (11.12) 

F F extension of total 

since A% = 1, a® = 0, atotai 7^ 0. (II. 7) is now proven. 
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Let us now check the other monomials in the v : s. If F = (Zi, . . . , l T ) is a forest we 
say that it creates the partition V of I n if V is the set of clusters of F. We then write, with a 
slight abuse of notation ax> = Yl{ij}/v u {ij}^ = Yl{ij}/v v {ij}' wnere {u}/^ means that i 
and j are in the same element or component of T>. 

Monomials generated by (II. 5) are of the form bx>, V being created by a sub-o- forest 
of F. Remark the if V is created by a forest of length r then \V\ = n — r. There are two 
cases to be treated: 
Case 1: V ^ {I n } 

Here the coefficient of bj) is zero in the left hand side, we must show that also for 
the right hand side. Let v = n — r, then (st abbreviates "such that") 



E 



— n V « - a£) 



~>V F=(l 1} ...,l T ) L L 

o — forest st 
(Zl , . . ■ creates X> 



o-forest st F=(F 1 ,F 2 ) is 

creating T> an _ forest 



Again, summation is over all possible r's. We now arrive at the heart of the inductive 
argument. We show using (II. 7) that F\ being fixed the sum on F 2 vanishes. In fact, this 
sum is the analog of (II. 7) when instead of I n we use V as a point set. If a and b are two 
elements of the partition V, we let U{ a b} = J2iea jeb u {ij}- Given an o-forest F = (F 1? F 2 ) on 
I n such that Fi creates V, F 2 = {l v +\, ■ ■ ■ ,l T ) induces an o-forest F = (Ii, . . . ,l T - u ) on V in 
the following way: if l K = {ij}, v + 1 < k < /i, with i G a and j E b, a and 6 elements of X>, 
then set l K - v = {ab}. F is simply obtained by forgetting the details of structure inside the 
components of V, which are due to the forest F\. Demanding that the whole of (Fi, F 2 ) be 
an o-forest insures that, in the process, F remains a forest. As a result, the analog of the 
for the forest F are = a^ +u — ax>, for 1 < p < fi — v\ and with f = r — v we have 



E (-D T - n^lfe) = E M) T ^ • Cn.14) 

^=(«„+i,...,«r) IW+ilS f=(7 1; ...J-) Ci- 



st F=(F 1 ,F 2 ) 
an o — forest 



The last sum is zero by (II. 7) because \D\ > 2, since we are in the case V ^ {I n }- 
Case 2: V = {I n } 



Here a-p = a tota i- Equating the coefficients of bx> = Yliev n Vl on both sides needs 

showing 

1= E 7 V?"' 1 ^ (IL15) 

F=(i 1 _« n _ 1 ) ( atotal -<)••• ( a totai - a r n _ 2 ) 

o — forest 

This last identity is shown by induction on n. Remark that here F is a complete tree covering 
I n . If n = 1, the only forest is the empty one, its contribution is the empty product i.e. 1. 

The induction step form n to n + 1 needs an argument similar to the former case: 

ui t . . . ui n _ 

F=(i7^., ln ) ( fltotal - • • • ( a total - fln-l) 

o — forest on I n 

Fi^fJi) atotal F 2 = (l7^.,i„) (°total - af ) • • • (atotai - a>n-i) 

1 1 st F=(F 1: F 2 ) is 



an o — forest on I 



n+1 



The link h = being fixed, this determines a partition J n = U {{i}/i G 

J n+ i, i ^ a, j ^} for which we can repeat the treatment of the partition V in the 
preceding case: introduce variables U{ a b} = J2iea jeb u {ij} f° r a an< l ^ n an< l induced o- 
forests F = (I l7 . . . , I n _i) on J n for any o-forest F 2 = (h, ■ ■ ■ , l n ) on sucn that ^ = (-^1, .P2) 
be an o-forest on I n +i- Remark that 

def 



tttotal 



= f E U Wb} = ( E ^ J ~ = fl total - af , (H- 1 ?) 



and for < p < n - 2, = aj +1 - af , therefore a to tai ~ a p = (atotai - af ) - (aj +1 - ai) = 



F 

atotai — °p+l- 



Hence 

uu . . . Ul 



E 



F 2=( i 2 ,. ..,,„) ( fl totai - ^ ) ■ ■ ■ (atotai - a n _i) 

st F = (F 1 ,F 2 ) is 



an o — forest on J 



n + 1 



Elty ... My 
J ^ (11.18) 

= (Ti,..J„_i) Atotai - a Q J . . . Atotai - a n _ 2 J 



F 

o — forest on J 



but the last sum is 1 by the induction hypothesis since \J n \ = n. Finally in (11.15) there 
remains 

= atotai =1 (ILig) 

, atotai atotai 
tl 

This completes the proof of Lemma II. 2 . ■ 



Proof of theorem II. 1: 

We can return now to the Brydges-Kennedy forest formula and first prove it for any 
complex numbers ui, I e V n such that for all subset X of V n , Xwex u i 7^ 0- In fact we can 
rewrite the right hand side of (II. 1) using o-forests instead of u- forests. Let $ be a u- forest 
and choose an ordering F = (Zi, . . . ,i T ) of its elements. We can slice the integral according 
to the ordering of the his thereby giving a contribution for $ 

E / fn^)fn^)exp(E^f(h).^) , (IL20) 

*e&r Jl>hl «(i) > - >hl «M >0 V=i / V=i / ier n 

& T is the permutation group on r elements. For any of its elements a let us denote the 
o-forest (Zcr(i), • • • , l a ( T )) by ^o-- Let us define, if F is an o-forest, the function h^y = h\ v 
where v is the largest index in the ordering provided by F of the links appearing in the unique 
path of F connecting % and j. If no such path exists we again put by convention h^-y = 0. 
Thus the contribution of 5 becomes 

E / (n^)(n^)exp(E^(h).^) (11.21) 

aeeT Ji>h^ (1) >...>h^ (T) >o\ u=1 / \ v=1 J l€Vn 

and the right hand side of (II. 1) 

1= E / (ri^ „)(n«Oexp(5>f(b).«0 • (n-22) 

f=(i 1 ,...,! t ) Jl>hi 1 >...>hi T >0 \ u=1 J \ v=1 J lePn 

o — forest 

Now we perform the change of variables t\ = hi x — hi 2 , . . . , t T -± = hi T _ 1 — hi T , t T = 

hi T then 

1= E u h ■ ■ ■ u i T t >0 t >0 dt i ■■■dt T exp(E^°v) • (11.23) 

F=(! 1 ,...,i T ) J t{+.'.'.'+tr<l V=\ 



jrest 



There shows up the when collecting in the exponential the it{y}'s multiplied by a given 
t v . Since s v = t v + t u +i + . . . + t T , the corresponding pairs {ij} are those for which all the 
links appearing in the unique path of F connecting i and j have an index at most equal to 
u, i.e. {ij}/u. Since ciq = 0, we can rewrite (11.23) as 



/= ui 1 ...ui T dt . . .dt T expfS^^af ) , 

^— ' I t >0,...,t T >0 ^^-^ " ' 

F=(! 1 ,...,i T ) ^ t + ...+t T =l I/=0 

o — forest 



(11.24) 



then from lemma II. 1 we obtain 



F 



/= y, "'■■■■^■(E n %.«-<) ' ' (IL25) 



F=(h,-M V=0 LL ^ V v v 

o — forest 



Then if we substitute for the indeterminates ui the complex numbers we have now 
(this is allowed since the factors in the denominators — are of the form ± J2iex u i ^ 0) 
and for the vi the complex numbers exp(itj), the lemma II. 2 just proves (II. 1). We can remove 
the condition Yliex ui ^ of being in a finite number of hyperplanes of (C Pn , by density and 
continuity of both sides of (II. 1). To prove it for an arbitrary commutative Banach algebra 
B we need only note that the identity we are about to show, is of the form f(ui, . . . , Uk) = 0, 
where the function / is analytic, its value for w's in a commutative Banach algebra being 
defined by its power series expansion. Thus proving it for all complex numbers is enough to 
entail the vanishing of the coefficients of the power series and to prove the identity for any 
u's in B. ■ 



Proof of theorem II. 2: 

We recall the notion of connected set function. If X is a finite set and ip is a map 
from V(X), the power set of X, to a commutative algebra £>, the connected function of ip is 
ip c : V(X) — > B defined inductively by 

k 

VCcX, V(C)= HMCi) , V> c (0) = i • (H.26) 

n={c 1 ,...,c fc } i=i 

partition of C 

Let X = I n , B be the field of complex numbers, V/ G V n ,ei=e Ul — 1, and let ip(C) = 
Yl{ij}cc(l + e {ij})- We ma ke also the assumption that for any non empty subset Q of V n , 
J2ieQ u i ® • K is well known [Bl] that 

mc)= E (IN ( IL27 ) 

25 graph connecting C lE& 

("connecting C" means that for every i ^ j in C there exists a path in going from i to j, 
the empty graph connects C if |C| < 1). 
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We claim that 

MC)= £ (n / ^)(n^)ex P ( J2 w mW- u m) ■ ( IL28 ) 

T tree connecting C^lEl '^lel ' ^{ij}CC ' 

This is trivial if \C\ < 1. If \C\ > 2, let T be a non trivial tree connecting C. We associate 
the sequence C % = (Cj)„ 6 ]N of disjoint subsets of C such that is the z/-th layer of X. 
Note that = {r c }, U^j^Cj = C, and 3^ G IN such that v < v% Cj ^ and 
z/ > zaj =>- CjJ = 0; a sequence C = (Cu) ue j<\ sharing those three properties will be called 
admissible. 

We will compute in the right hand side of (11.28) the contribution of the trees 
corresponding to a given sequence C % = C. Let {ij} be a pair in C with i G C^, j G C„. and 
i^j > Uj. The only way ty^-^(w) can be non zero is that Vi = Uj or Vi = Uj + 1. For z/^ = za,- 
we just get a factor e u ^> = 1 + since w^}( w ) = 1- 

But if Vi = Vj + 1, let i' be the element of C v . such that G X; is the ancestor 
of z. Let us perform the integration explicitly in The pairs / such that wf(Mv) d = w^niy 

are the pairs {ij} with j G C„ { _i, the corresponding contribution is 



/ M{ tf }.exp > u {y} =u {ii , } . = . (II. 

Jo V • eCi/ ._ i / ^iea.-i^ii} 



29) 



Now, building a tree with prescribed C T means linking the points of Cf to the root 
rc, then for each i in choosing an ancestor z' in , then for each z in Cjf choosing an 
ancestor i' in Cj, and so on, until we exhaust C. The choice of ancestor for points in the 
same layer are completely independent operations. Summing the contributions (11.29) for 
all choices of ancestors for i is simply 

u i^'} ■ v ^ = exp l 2^ J ~~ 1 



= n ( i +^})- i = e (n^o • ( il3 °) 

Therefore the right hand side of (11.28) is 

e n((n( e n^i))x( n a+q«))) . aum 



admis 

sequence C ,7^0 
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and is to be compared with (11.27). 

Note that, in a connected graph (J5, we can define the distance d<g,{i,j) between two 
points i and j by the minimal number of links in a path in (5 going from i to j. Given 
a preferred point in C, namely the root rc, and a graph (25 connecting C, we can define 
an admissible sequence C® = (C®) ue j$, by letting C® d ={i G C/d&(i,rc) = v}- This is 
consistent with our previous definition since a tree is obviously a special case of graph. 

We just have now to convince the reader that for a given admissible C 

n((n( e n«<«))*( n a +««»>)) = e it< o^* 

connecting C 
such that C®=C 

Note that building a graph (3 with = C is first linking the points i of Ci to the root rc, 
then for each % in C2 choosing a non empty subset Ji C Ci of points j to which z will be 
linked, then for i in C3 we again choose a non empty subset C C2 of points j to be linked 
with, and so on. Having built such a "skeleton" we, as a final step, have complete freedom 
to add all the links {ij}, for i and j in the same C v , that we want. Remark that in such 
a graph (5, there is no link {ij} with % G C Ui , j G C Vj and Vj + 1 < Vi, for we would have 
d&( r c, i) < d®(rc,j) + 1 = ; + 1 < Vi = d®(r c , i), a contradiction. 
Now (11.28) is proven, so the right hand side of (II. 2) is 

k 

n={d,...,c fc } i=i ; 6 p„ z 6 P„ 

partition of 

Finally we can remove the conditions XlzeQ u z an d take the ui in any commutative 
Banach algebra £>, as we did for the proof of theorem II. 1. ■ 

III. A first Generalization: The Taylor forest formulas 

In practice we not only need "algebraic" cluster expansions, but also "Taylor" cluster 
expansion, in which one typically makes a Taylor expansion with integral remainder to inter- 
polate between a coupled situation (parameters set to 1) and an uncoupled one (parameters 
set to 0). 
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For instance a "pair of cubes" cluster expansion [Rl] boils down to applying the 
following interpolation formula: 

where /(a?i, . . . , x n ) is a smooth function on IR n , and is the vector of coordinates 

(xi, . . . , x n ) defined by Xi = if % I and Xi = hi if % G I. Such a formula is easy to prove. 

The "Taylor" tree cluster expansion describes a similar interpolation formula. Let 
S be the space of smooth functions from IR 77 ™ to an arbitrary Banach space V. An element 
of IR 77 " will be generally denoted by x = (xi) leV . The vector with all entries equal to 1 
will be denoted by 1. Applied to an element H of 5, we can state two different Taylor forest 
formulas depending on which of theorem II. 1 or II. 2 we use for its derivation. 

Theorem III.l (The Brydges-Kennedy Taylor forest formula) 

H W = £ (nMftnlrHWM) • (m.2) 

$ u-forest \es y v y 

f/ere X~ K (h) is the vector {xi) leV of IR 77 ™ defined by x\ = hf(h), which is the value at 
which we evaluate the complicated derivative of H. 



and 



Theorem III. 2 (The rooted Taylor forest formula) 

h w = e (n /' dw ) ((n 4r) H ) (w) . (111.3) 

5 u-forest \es y x ies 7 

ffere X£(w) zs £/ie vector (xi) leVn ofiR Vn defined by x\ = wf(w), which is the value at which 
we evaluate the complicated derivative of H. 

Proof: These formulas look more general than the algebraic forest formulas (II. 1) and (II. 2) 
(obtained by letting H(x) = exp(^ ;e7 , xiui) ) but are in fact a consequence of them. For 
any integer p > 0, let S p be the space of polynomial functions in the xi's of total degree at 
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most p. S p can be dressed in a finite dimensional Banach space structure. Let B p be the 
algebra of linear operators on S p generated by the derivations d/dxi for I e V n , equipped with 
the operator norm. B p is a finite dimensional commutative Banach algebra. If we consider 
the ui d = d/dxi in B p we have the right to use the multiple forest formula of the preceding 
section. 

We then remark that, in S p: ex.p(Xd/dxi) is the translation operator by the vector 
—Xei, where {e{) l£V is the canonical base of IR/^™ and the exponential is defined by its power 
series, which converges because d/dxi is a bounded operator on S p . It reduces to a translation 
simply because Taylor's formula is finite and exact for polynomials. 

Since exp{J2iev n Xid/dxi) is the translation by A = — J2iev n ^ iei i- e - the operator 
on S p which maps the polynomial function H to H' : x i— > if (x — A), we have proven (III. 2) for 
any H in S p . But p is arbitrary so we have proven the formula for any polynomial function H. 
By density of the polynomials in the Frechet space S for the C°° topology and the continuity 
relatively to that topology of both sides of (III. 2), we extend the validity of the equation to 
all H in S. We have proved formula (III. 2) but the argument is the same for formula (III. 3), 
just replace "h" by "w" . ■ 



IV. A second generalization: The jungle formulas 

Let m > 1 be an integer, and (uf) leV , 1 < k < m, be m families of n(n — l)/2 
elements of a commutative Banach algebra B. An m-jungle is a sequence T = (5i, • • • , dm) 
of u- forests on I n such that JiC.C 5m- 

The multiple forest formula states: 

Theorem IV. 1 (The algebraic Brydges-Kennedy jungle formula) 

ex P ( £ «f)= £ (n r*-)(n( n «?)) 

l<fe<m m— jungle 

• «-p(EE^( h )-«f) ^v- 1 ) 
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where $o = by convention, h is the vector (hi) le g m and the functions fr^V.j(h) are defined 
in the following manner: 

- If i and j are not connected by $k let ftjV j(h) = 0. 

- If i and j are connected by $k but not by Bk-i let 

h^f } (h) = M [h h l e L Sk {z 3 } n (3*\3*-i)} 

(recall that L${ij} is the unique path in the forest 5 connecting i to j). 

- If i and j are connected by $k-i let fr^V* (h) = 1. 

Proof: By induction. The case m = 1 was treated in Section II. For the induction step 
from m to m + 1, we sum over the last forest $ m +i' 

e (n Mffff n «f))«p(EE*r»-«f) 

(m + l) —jungle 

= e (iijM(n( n «.*)) «p(e e *rv).«f) 

^'=(Si, ..,5m) V ZeS m 70 7 \fc=l ^e^fcXSfe-l 7 / V fc=li6P„ 7 

m— jungle 

E ( II /^H II ^expf^^lh).^ 1 ) (IV.2) 



5m C 5 m |l 



where h = (^i)/ e y m+l) h' = (hi) le $ m and we have noted that if 1 < < m then (h) = 
' (h'). To perform the summation over 5" m +i, we will use the forest formula of section 
II and our favorite argument of forgetting the details of the tree structure up to $ m , to 
concentrate on what dm+i brings as new connections between the existing clusters. We 
introduce the partition V of I n created by $ m = {h, . . . , l„} and the u-forest on V, dm+i = 
{7i, . . . , l T - v } induced by 3m+i\3m = . . . ,l T } with v < r. The definitions are the 

same as in the proof of Lemma II. 2 except that we have u-forests instead of o-forests. 

For a link {ab} between two elements a and b of V, let U{ a b} = Yliea jeb • 
Summing over 5m+i, u-forest on I n containing # m , with the "propagators" u^Jy is the same 
as summing over the u-forests $ m +i on V with the "propagators" U{ a b} i-e. 

( II / <%)( II ^) eX p( J2 ^{abV^-Hab} 

15 



e ( n M( n «r +i W( e fcsr'w--^ 

such that f m Ci m |l Ujit'D 

(IV.3) 

where {ij}^D means "i and j are in different components of P". By the forest formula of 
section II, the left hand side of the last equality is 



exp Hab}) (IV.4) 



{ab}(ZT> 



that is 



ex pf E u m) ■ ( IV - 5 ) 



The right hand side is almost the partial sum we want to perform on $ m +i in (IV. 2) 
there misses 



H E *«r 1 (»)-«K 1 ■ < IV - 6 > 



. {ij}ev n 



The sum is over all pairs {ij} in I n such that {ij}/V i.e. i and j are connected by # m . But 
then the definition of the functions hf' (h) tells us ft^V™ +1 (h) = 1; so the missing factor 



becomes 



ex pf E u m) ■ ( IV - 7 ) 

\{ij}e7>„ / 



In conclusion 



E ( II f<"")( II «r +1 )ex P (^^"' +i (h). M r«) 

i-forcst on I„ M^'S . . \ ir " / M/^Tr . \ Tr / McU ' 



S m+ l u-forest on J„ \edm+l\dm U ' ^Sm+AZm ^Tn 

St 5 777, C "15 777, -|_ X 



= exp( ■ 



(IV.8) 



The sum over the m-jungle T' is now exp(5^^_ 1 Xwep uf) by the induction hypothesis, and 
the right hand side of (IV. 2) is finally summed to exp(Y^)™j~ 1 L J2 le -p n uf) as wanted. ■ 
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We can generalize in a similar way the rooted forest formula. Given an m-jungle 
J~ = (3l> • • • 7 3m), we introduce the notation w for the vector { w l)ie$ m i and iu^V.j(w) for the 
functions defined by: 

- If % and j are not connected by 3fc let w^jy ( w ) = 0. 

- If % and j are connected by 3&-1 let tu^Vj(w) = 1. 

- If i and j are connected by 3& but not by 3fc-i, let rc be the root of the cluster 
of 3fc containing z and j, define lk(i) to be the number of links of 3fc\3fc-i in the unique 
path that goes from % to rc, and similarly for j. Now if |/fc(i) — /fc(j)| > 2 put iu^V.j(w) = 0, 
if /fc(i) = / fc (j) put w^j(w) = 1, and if l k (i) = h{j) + 1 for instance, take w^(w) = wi 
where / is the first link of 3fc\3fc-i on the path that goes from i to rc- It is about the same 
definition as for theorem II. 2 but we take only into account the links in 3fc\3fc-i- 

We can now state 

Theorem IV. 2 (The algebraic rooted jungle formula) 



exp 



e «?) = e (n M(n( n «o 



le-Pn ' N»l »m) M€3m" 7 V fc = l JSSAfffc-l 

1 < fe < m m— jungle 



• ex p(E E w f ,k w- u i) ■ ( IV - 9 ) 

Proof: The deduction of the jungle formula from the forest formula is the same as in the 
Brydges-Kennedy case. Remark that the property that was used was the following. Once 3m 
is fixed as well as the partition V of I n it creates, the numbers w^™ +1 (w) are unchanged 
if we allow % to travel freely in the component eij of D it belongs to, and the same for j in 
cij. Furthermore w^™ +1 (w) = ^^^(w) where 3m+i is the forest on V induced by 3m+i, 
and w = (wjjj^ is defined by Wj = wi where / is the unique link in $ m +i\$m inducing 
the link I in V. The functions «7?(w) are defined by the same algorithm as in theorem II. 2 
but with V instead of I n as a point set, and with the following induced choice of root. If C 
is a non empty subset or cluster of V, C = {a±, . . . , a^}, let C = \J^ =x a v . It is a non empty 
subset of I n and already has a chosen root rc- Then 3! v, 1 < v < fx such that rc G a v . So 
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take r-£ = a v G C to be the induced root of C, that will get involved in the definition of the 
Wj(W). 

Finally, and this is the last ingredient that allows us to copy the precedent proof, 
we have forced w^ i 'J r ^ +1 (w) = 1 if i and j fall in the same component of V so that we once 
again recover the missing factor 




(IV.10) 



■ 

The reader must have noticed along this proof that we used a fairly general recipe 
to obtain a jungle formula from a forest formula. We also have a recipe to get Taylor 
interpolation formulas from algebraic ones: simply apply them to differentiation operators 
for the M/'s. 

Let V n ,m = V n x I m and S be the space of smooth functions from IR ™ m to an arbi- 
trary Banach space V. An element of IR^™" 1 will be generally denoted by x = (xf)^ k ^ eV 
The vector with all entries equal to 1 will be noted 1. H is an arbitrary element of S. We 
can state the following theorems, whose proofs are rewritings of that of theorem III.l. 

Theorem IV. 3 (The Brydges-Kennedy Taylor jungle formula) 

e (n M((n( n i*)) H )( x ° K ^ ■ ^ 

^=(*i,...,* m ) Hedm J0 7 V v fc=i iedk\dk-i 1 7 / 

m — jungle 

Here X^ K (h) is the vector (xf)^ lk ^ eV ofiR Vn ' m defined by xf = hf' k (h), which is the 
value at which we evaluate the complicated derivative of H. 

and 

Theorem IV. 4 (The rooted Taylor jungle formula) 
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Here X£(w) is the vector (x k )^ lk ^ eV of TR Pn ' m defined by xf = wf ,k (w), which is the 
value at which we evaluate the complicated derivative of H. 

Let us conclude this section by recalling the important positivity property of the 
Brydges-Kennedy forest and jungle formulas. 

Theorem IV. 5 (Positivity) Let T be a m-jungle, and M k , 1 < k < m be any sequence 
of m positive symmetric n by n matrices with entries m k j . Then the interpolated matrices 
M k (h) with entries m^-(h) = hfj ,h (h)m k j if i < j and m^(h) = ra\ are all positive. 

Proof: The definition of M k (h) only involves the forests $k-i and $k and the h^s for / in 
Note that, for every link /, hf ,k (h) = hf k (h), the right hand side is the hf k (h) function 
of the simple forest formalism of Section II, where h = (hi) l€ ^ with hi = hi if I e $k\-5k-i 
and hi = 1 if I e $k-i- In fact, for all points i and j connected by $k-i, i-e. such that 
L$ k {ij} C fffc-i, we have 

def 

hf"(h) = mf{h h I G L dk {ij}} = inf{l} = 1 = hf' h (h) . (IV.13) 

So we just need to prove the positivity in the simple forest formalism. This has 
already been done by Brydges in [B2], but we give here a proof along the lines of our derivation 
of the formula (II. 1). 

Suppose we have a forest J with parameters h = (hi) le ~ between and 1, and 
a positive symmetric matrix M with entries m^- that is interpolated by M(h) with entries 
Tfiij (h) = mijhf(h). Let us number the links of $ into an o-forest F = (h, . . . .l T ) in order 
that 1 > hi x > . . . > hi T > . We perform the same change of variables as in the formula 
(11.23). If i ^ j let us denote by U^jy the matrix with zero entries except at the locations 
and (j, i) where we put 1. Let Mdi ag be the matrix with entries maSij the diagonal part 
of M. Then 

M(h) = M diag + hf(h).Uimi (IV.14) 

iev n 

and for the same reason as (11.23) 

T 

M(h) = M diag + ( IV - 15 ) 
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where = u {ij} m ij, and finally 

T T 

M(h) = (l-J2tv)M diag + J2 t AM? + M diag ) . (IV.16) 

Now note that for 1 < v < r, t v > and X^=i ^ v = h\ x < 1 so 1 — 5X=i ^ ^ 0' 
and Mjf + Mdi ag is a positive symmetric matrix. Indeed let us fix v and consider the o- 
forest F v = (Zi, . . . , l v ) completed at stage v from the given o- forest F. We can reorder the 
base vectors so that M looks like a matrix with blocks corresponding to indexes in the same 
component of F v . Taking only the diagonal blocks and putting the others to we get exactly 
the matrix + Mdi ag . It is clear that such a process preserves the positivity of the matrix. 
M(h), a sum of positive matrices with non negative multipliers, is positive. ■ 

Note that the rooted formulas do not preserve positivity as can be seen on simple 
examples, with 4 points for instance. 

V. Some concrete examples 

V.A) Gaussian measures perturbed by a small interaction 

In constructive field theory cluster expansions are typically used to perform the 
thermodynamic limit of a field theory with cutoffs and a small local interaction Ia(4>) i n a 
finite volume A. The cluster expansion expresses the partition function as a polymer gas with 
hard core conditions (see e.g. [R2]). The set I n is then made of a partition of A into (unit 
size) cubes, and clusters are subsets of such cubes. 

For instance consider the free bosonic Gaussian measure d\xc i n IR d defined by a 
covariance C with ultraviolet cutoff and good decrease at infinity. A standard example is 

so that 

( P 2 +m 2 ) 

C(p) = — 9 =- (V.A.2) 

We perturb this free theory by adding an interaction such as e 9 §^ ^dx, and we 
want to perform the thermodynamic limit A — > oo, that is to define and study an intensive 
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quantity such as the pressure 

p= lim i-logZ(A) , (V.A.3) 
A^IR l A l 

where the partition function Z(A) in a finite volume A is 

Z(A) = J df, c (<P)e- g ^ 4,4{x)dx , (V.A.4) 

Let us explain how the Taylor formula (III. 2) performs the task of rewriting the 
partition function as a dilute gas of clusters with hard core interaction. We write A = 
Ui e i n bi, where each bi is a unit cube, and define Xb as the characteristic function of 6, and 
Xa = J2iei n Xh- Since the interaction lies entirely within A, the covariance C in (V.A.4) can 
be replaced by Ca = xa(x)C(x, y)xA(y) without changing the value of Z(A). Moreover Ca 
can be interpolated, defining for / = {z, j} e V n 

n 

C A ({xi) lGPn )(x,y) = J2xb i {x)C(x,y)xb i (v) 

i=i 

+ Yl x {ij}{Xb i (x)C(x,y)xb j (y) +Xb j (x)C(x,y)xb i (y)) (V.A.5) 
{i,j}er n 

Remark that Ca(1,...,1) = Ca- Now we apply the Taylor formula (III. 2) with the func- 
tion H being the partition function obtained by replacing in (V.A.3-4) the covariance C 
by CA((xi) leVn ). Here it is crucial to use the positivity theorem IV. 5, in order for the in- 
terpolated covariance to remain positive, hence for the corresponding normalized Gaussian 
measure to remain well defined. From the rules of Gaussian integration of polynomials, we 
can compute the effect of deriving with respect to a given xi parameter, and we obtain that 
(III. 2) in this case takes the form 

Z(A) = H(l) = W c^ CA(xr(h)) (0)(n f dhi) 

{ || /4 t «»,w» J (»)c(«)^4K sJ *''" >i * <V ' A ' 6) 

Since both the local interaction and the covariance as a matrix factor ize over the clusters of 
the forest the corresponding contributions in (V.A.6) themselves factorize, which means 
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that (V.A.6) can also be rewritten as a gas of non-overlapping clusters, each of which has an 
amplitude given by a tree formula: 

/n 
d^ CA (<P)e lA ^ = £ RAW (V.A.7) 

sots {Yi,...,V„} i=i 
Y i nY J =0,UY i =A 

A ( y )= e (n / / ^ ( xi-(h))(0) 

T on y Z6l^° ^ 

{ n /^(.te.W,*)^).-^* (V.A-8) 

2 — { } £ X 

where 6j and are the two ends of the line /, and the sum is over trees X which connect 
together the set Y, hence have exactly \Y\ — 1 elements (if \Y\ = 1, X = connects V). The 
measure oI\i Cy ^ x bk (h)})(0) is the normalized Gaussian measure with (positive) covariance 

Cy(X| x (h))(x, 2 /) = xy(x)(^(h)(x, 2 /))C(x,y)xy(y) (V.A.9) 

where /i<j(h)(£, y) is 1 if x and y belong to the same cube, and otherwise it is the infimum of 
the parameters hi for I in the unique path L<x(b(x), b(y)) which in the tree X joins the cube 
b(x) containing x to the cube b(y) containing y. 

Formula (V.A.8) is somewhat shorter than the different formulas of [Rl-2], and can 
be used in the same way to check that given any constant K 1 for small enough g with Re g > 

\A{Y)\R\ Y \ < 1 (V.A.10) 

Y st OeY 

Proof of (V.A.10) requires the slightly cumbersome computation of the action of the functional 
derivatives in (V.A.8) and a bound on the resulting functional integral. The method is 
identical to [Rl-2]. Remark that although the full amplitudes A(Y) defined in (V.A.8) must 
be identical to those in [Rl-2], the subcontributions associated to particular trees are different. 

The Mayer expansion defined below allows to deduce from (V.A.10) the existence 
and e.g. the Borel summability in g of thermodynamic functions such as the pressure p 
defined by (V.A.3). 

B) The Mayer expansion 
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In the cluster expansion (V.A.7), the condition that the disjoint union of all clusters 
is A is a global annoying constraint. Remark that the polymer amplitudes are translation 
invariant. In particular the trivial amplitude of a singleton cluster Y = {b} is a number Aq 
independent of b. Redefining A r (Y) = A(Y)/A 1 ^ 1 and Z r (A) = Z(A)/A^ 1 we quotient out 
all the trivial clusters so that 

n 

z r (A)=i+Yi e n^rco ( v - b - 1 ) 

n >l sots {Y 1 ,...,Y n } i=l 
\Y t \>2 , Y-nY^fl 

This is the partition function of a polymer gas: the sums over individual polymers would 
be independent were it not for the hard core constraints Yj fl Yj = 0. Adding in an infinite 
number of vanishing terms, we can replace the sum in (V.B.I) by a sum over ordered sequences 
(Yi, ...,Y n ) of polymers with hard core interaction and a symmetrizing factor l/n\ coming 
from the replacement of sets by sequences. 

^( a ) =i +£^ e fiMYi) n e ~ v{Yi ' Y]) ( v - b - 2 ) 

■ -I 

n>l sequences (Yi,...,Y n ) t=± l<i<J<n 

where the hard core interaction is V(X, Y) = if X fl Y = 0, and V(X, Y) = +oo if 
X fl Y 7^ To factorize again this formula we cannot apply directly the algebraic Brydges- 
Kennedy formula (II. 1), because the interaction V can be infinite, but the Taylor forest 
formula (III.2) easily does the job. More precisely we define now I n as our set of indices, 
and define e^.j = [e~ v< ^ XuX ^ — 1), for i ^ j. For a fixed sequence (Yi, . . . , Y n ) of polymers, 
consider the function 

H((x l ) leP J= Y[(l + Xl eY) (V.B.3) 
iev n 

Rewrite (V.B.2) as 

1 

Z r (A) = l + ]T- H(i)HMK) 

n>l sequences (Y 1 ,...,Y n ) i = l 

\Yi\>2 

n >l ' sequences (Yj , . . . , Y„ ) i= 1 § J ° ' \ 6 £ ' ' 

\Yi\>2 
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= E^(Eh E (n^«))c , (*i,...,n))" (v.b.4) 

n>0 fc>l sequences (Y± , . . . ,Y k ) i = \ 

mi>2 

where 

c T (Y U ...,Y k) = y: ik 

G connected graph 
on{l,....,fc} 

e (n/ 1 n e n (^iw) (v.b. 5) 

X tree on {l,....,fe} ^ZeX ^° ' NePfe, «€l ' ZeP fc , 

where /if (h) is, if / = {i?}, the infimum of the parameters hy for V in the unique path L%{ij} 
which in the tree X joins % to j. 

We obtain immediately that 

logZ r (A) = ^l £ rn^(^))c T (^i,-,n) (V.B.6) 

k>l sequences (Y± , . . . ,Y k ) \ = 1 ' 

im>2 

Formulas (V.B.5-6) are more explicit than those used in [Rl-2] and have all desired advantages 
(every tree coefficient forces the necessary links and is bounded by 1, since |(1 + hf (h)e^)| < 
1). They can therefore be used together with (V.A.10) to control in a similar way the 
thermodynamic function p = Aq + lim A ^-jr^d logZ r (A). 

C) A single formula for the succession of a cluster and a Mayer expansion 

Formula (V.B.6) involves trees and forests on two different kinds of objects. First 
are the links between boxes in our lattice that make the clusters YJ. Then we have Mayer 
links between indexes i in I k . This leads to complications, if we iterate the process, when 
doing a renormalization group analysis. So we propose a formula with two types of links 
but on the same objects, it is a 2-jungle formula. For simplicity we write it in the algebraic 
setting of section II, i.e. with Z = exp(^-uj) instead of Z r (A), and perform the apparently 
stupid operation logexp(^-u;) ! 

The links / are between elements of 7 n , but now we introduce a new set of objects 
V n = I n x IN. If I is a link {(&, k), (b', k')} between different elements of T) n we let uj — U{bb'} 
and 5j = dbb' where the second delta is Kroneker's. Let T = #2) be a 2-jungle on T> n , we 
say that it is nice if it fulfills the following requirements: 
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- 3i is only made of horizontal links i.e. of the form {(&, k), (6', k)}; 

- $2\5i is only made of vertical links of length 1 i.e. of the form {(&, k), (6, k + 1)}; 

- #2 is a non empty tree; 

- there are no singletons among the clusters of 5i in the support of #2; 

- there is a unique cluster of #1 in the support of #2 that lies in the 0-th level 



- if we fix a root r of #2 in that unique cluster, then for every element (6, k) in 
the support of $2, k is the number of vertical links on the path going from (6, k) to r (this 
property is independent of the choice of r). 

We denote by k{!F) the number of clusters of #1 in the support of #2- We introduce 
the following lexicographical order on V n , (6, k) -< (&', k') if and only if k < k' or (k = k' 
and b < b'). This is the order we use to choose a root for each non empty finite subset of 
V n according to the rooted jungle formalism. We let h = (hj)j e ^ i and w = ( w i)j e g 2 \g 1 an d 
we take the compound (h, w) instead of the h or the w vectors that were used in (IV. 1) and 
(IV. 9). The functions hj ' and Wj ' are however defined exactly in the same manner as in 
section IV. We now claim that 

Theorem V.C.I 



/„ x {0}; 



logZ = 





(V.C.I) 



Proof: 



For m G IN, let V 



I n x {0, 1, . . . , m} C V n . As in (V.B.6) we have 



log Z = lim L 



rn 



(V.C.2) 



where 




(V.C.3) 
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and 

A{ y)= E (n«' /^WfE h f( h ) u i) ■ ( v - c - 4 ) 

X tree on Y Mel ^° ' M in Y ' 

This time we compute C T (Yi, . . . , Y 8 ) thanks to the rooted forest formalism. Note 

that 

n ( i+e k)= n n a-^w 

{ij}Cl n {ij}<Zl n {b,b')eY i y.Y j 

with all the x^jy set to 1. Then we apply the rooted Taylor forest formula III.3 and collect 
the connected parts. We get 

C T (Y 1 ,...,Y S )= J2 C(s;Y 1 ,...,Y s ;1) (V.C.6) 



% tree on /„ 



where 



c( s; yi,...,y s ;X)= (u [ dwi 

\ez Jo 

n ( e (-^) n (i-^}(w)5 CC o) 



(c,c')#(b,6') 



• n ( n (i-^cw)^)) • (v.c.7) 

The only trees X giving a non zero contribution are those where for any i ^ j in 
the same layer of X, Y~j fl Y} = 0. In fact for such i and j, {ij} X and w fijy( w ) = 1 so in 
(V.C.7) we get a factor 

{l-8 w ) = e- v <?*™ (V.C.8) 

(b,b')eY,xYj 

forcing the non overlapping condition. 

Given a sequence (Yi, . . . , Y s ) and a tree X satisfying that property, we can define 
for i e Ik, Yi = Yi x C V n ^ m . Here is the height of vertex i in the tree X. The 

Yi are all disjoint. 

For a finite subset F of V n let 

3(F) = £ exp(_E/f (h)« 7 ) (V.C.9) 



X tree on Y lei l&Y 
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so that A(Yi) = A(Yi) for every i e I s . Now define the (unordered) set O = 
0(Y 1 , . . . ,Y S ; X) = {Yi,...,Y s }. It is easy to see that O is a non empty set of disjoint 
polymers (|Y| > 2 for each Y G O) lying in T> n ^ m , with exactly one element in the ground 
level I n x {0}, furthermore, the labels k of the occupied levels I n x {k} form an interval 
{0,1, . . . ,q}, q < \0\ — 1. An (9 verifying these properties is called admissible. 
We will sum over admissible O's in (V.C.3) 

£m= E |^r(n^ F )) E c^y!,...,^;^ • (V.C.10) 

admissible O ' Vpfl (n. ■■-.*«);* 

st o=o(y 1 ,...,y s ;'r) 

Knowing the sequence Yi,...,Y s ), 1 naturally induces a tree on (9, X = 
X(Yi, . . . , Y s ; X) by the rule {ij} E X if and only if {Y t ,Yj} e X. The tree X has a nat- 
ural root i.e. the unique Y of O lying in the 0-th level, moreover, for each Y the label of the 
level it belongs to is just its height in the tree X with the mentioned choice of root. Such a 
tree will also be called admissible. 

Let for any tree X on O 

u(o,%)=i[J\w L (il ( e nji-in^^)) 

Lei ° ^{Y Y'}£% (b,b')eYxY' (c,-')ey_xj' 

(c,c')^(i>,!>') 

• n n ji-in^r}) • ( v - c - n ) 

and 

iV(£>,X) = E 1 • (V.C.12) 

(Y 1 ,...,Y S );% 
st = 0(Yi,...,y s ;T) 
and % = %(Y 1 ,...,Y s ; l Z) 

It is clear now that 

L m = E |^t(II^( F ))_E U{0,%).N{0,%) . (V.C.13) 

admissible O ' VeO Ttree on 

Let us admit for the moment the following 

Lemma V.C.I The combinatoric factor N(0, X) is (\0\ — 1)1 for any admissible O and X. 
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Then 



L m = E W\(U A ^) E • (V.C.14) 



admissible o VeO X tree on 

T tree on O 



Besides, 



n = e ( n u j I dh i) ex p (e h i 1 ( h )^) ' ( v - c - is ) 

where we sum over all forests 3i made of horizontal links uj such that the connected compo- 
nents that are not singletons are exactly the elements of O. 

Finally remark that summing over trees X on O and, for every {Y, y'} e X, over 
(6, b ) e Y x F is the same as summing over forests $2\5i °f vertical links — 5j connecting in 
a tree the clusters Y £ O formed by It is now a matter of (rather unpleasant) checking 
the definitions to convince oneself that (V.C.14) is just the result of formula (V.C.I), except 
that the clusters are confined to X> n , m . However taking the limit m — > +oo removes this 
restriction. ■ 



Proof of lemma V.C.I : 

Since O and X are admissible, there exists a sequence Yi, . . . , Y 8 and a tree X such 
that O = 0(Yi, . . . , Y s ; X) and X = X(Yi, . . . , Y s ; X). In fact suppose we are given a one-to- 
one map a : O i— > I s such that cr(yi ) = 1, where Vi is the unique element of O lying in the 
ground level of V n . We can construct a sequence (Y^ 7 , . . . , Y S CT ) and a tree X' 7 fulfilling the 
previous requirements. If 7r denotes the projection of V n on J n , we let = n (Y) for all 

FgO, and X°" be the tree on I s induced by X through the correspondence a. 

It is easy to see that every pair ((Yi, . . . , Y s ); X) we are counting is obtained that 
way. The interesting thing is that there is a unique a giving this pair. If a and r are two 
maps such that ((Yi CT , . . . , Y s a ); X CT ) = ((Yi T , . . . , Y S T ); X r ), cr o r -1 must preserve the tree 
X°" = X T , but a o t -1 (1) = 1 so it preserves also the root, as a consequence it leaves the 
layers invariant. But, the corresponding Yj are disjoint contained in the same level of V n , 
thus the projections Yj are distinct and Y aOT -i^ = Yi forces a o r _1 (z) = % inside each layer. 
In conclusion a o r _1 = /(i. N(0, X) is the number of cr's that is (s — 1)!. ■ 
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If after all that the reader is somewhat sceptical about formula (V.C.I), he might 
find it amusing to check it for n = 2. 

Now to use formula (V.C.I) for more clever applications than computing the trivial 
expression logexp(^ m), we state the analogous result for logZ(A) in theory. We intro- 
duce a copy of the field (fi k for each k G IN, and write $ for the system (<^>fc ) fce j\[ - Then we 
define the interaction in S 1 ^), the support of $2, as 

e i(S(d2))(fy = Yl e - x fb* t{x)dx (V.C.16) 
(b,k)eS(S2) 

where we have identified the boxes b in A with their labels in I n . To a link I = {(&, k), (&', k')} 
we associate the operator 

Cj = 6 kk , fdxfdy C(x,y) 6 6 (V.C.17) 
Jb Jv O0k{x) d(p k >{x) 

that will play the role of uj. 5j is defined in the same fashion as before. Given a 2-jungle 
J= = (# l7 5 2 ) denote by S k ($ 2 ) the fc-th slice S($ 2 ) H (I n x {A;}) of 5(ff 2 ). We define the 
following Gaussian measure on the fields 4>k\s k {$ 2 ) suc h that k is an occupied level of #2 (i.e. 

S fc (S 2 )^0) 

^' h ' w (%(* 2 )) = II ^S' h,w '*(0*|s fc (ff a) ) (V.C18) 

occupied k 

where dfj,Q' h ' vv ' k ((j) k \ Sk (g 2 - ) ) has covariance C :F,h,w,k (x 1 y) = C(x, y) if x and y fall in the same 
b with (6, fc) G S k ($ 2 ); and C7^ h - W ' fc (x,y) = C(x, y).h^ kh{b , k , )} (h, w) if x G b, y G 6', 6 ^ 6' 

(M)eSk(?2) and (&',fc) 

We can now state a formula which performs a cluster and a Mayer expansion in a 
single move (recall that |A| = n) 



Theorem V.C.2 

-|S(5 2 )| / /•! \ / rl 



iogz(A)=niogA + ^ -w n / n (-#/ ^ 

nice 2-junglc lfctf2\tfl 

J] (l-f^wl^l^^^^lln^e 1 ^'' • (V.C.19) 



Proof Use first log Z(A) = nlogA + logZ r (A), then repeat exactly the same line of argu- 
ments than for Theorem V.C.I. ■ 
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From Theorem V.C.2 one can construct a series for the pressure limA^oo |A| _1 log Z(A) which 
is the sum of the trivial term logA plus a sum over nice 2 jungles that extend horizontally 
over the infinite lattice of all cubes covering IR d , and such that the well defined unique root 
of $2 is the particular cube containing the origin. (This requires to neglect boundary terms 
of order |A|( d_1 )/ d in logZ(A)). For small A this series is absolutely convergent. 

D) The resolvent expansion 

In many situations physical situations (polymers, disordered systems) we want to 
compute not a full theory but a single Green's function which is expressed as the inverse 
of some operator. This mathematical situation is formally equivalent to 0-component field 
theory (or to supersymmetric theories) in which the usual normalizing fermionic or bosonic 
determinants have been cancelled out. 

Consider a finite dimensional operator M on lR n with matrix elements rriij, and a 
norm strictly smaller than 1. The operator A = 1 _^ M is well defined through its power series 
expansion. Define M(x) as the matrix with elements x^yrriij if i ^ j and diagonal elements 
ma equal to those of M, and M^jy as the matrix with elements except (M^y) = rriij and 
(Mfyjy) .. = rriji. The Taylor forest formula (III. 2) applied to the operator H(x) = 
gives 

a = ttm = = e (n f ""•) (n / ) 1+A f(x S (h» +Eie8 x,M, 

(V.D.I) 

where the action of multiple derivatives is for compactness rewritten as a multiple Cauchy 
integral, and the analyticity radii in (V.D.I) are small enough. This does not look like a very 
clever rewriting, but it allows to reblock the power series for y+m accor ding to the set of 
sites truly visited in this power series. In particular if we impose a given entry to the matrix 
A the forest formula (V.D.I) reduces to a tree formula: 

^ = (t^) - E E tilt*'* 



1 + M J ^ ^ V Jo 

'13 yc{i,...,n) X tree on Y 
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where ly is the identity on R Y , My(Xj(h)) is the matrix on R Y with diagonal entries ma 
(for z G y) and off diagonal entries ^(h)mjj (for j £ 7, j £ Y), and My^ is the restriction 
of Mi to indexes in Y . 

This kind of formulas, eventually in combination with a large field/small field anal- 
ysis should be useful for interacting random walks [IM] and presumably for the study of 
disordered systems. 

E) Cluster expansions with large/small field conditions: the m = 2 jungle. 

Jungle formulas with several levels are interesting when there are various type of 
links with priority rules between them. The simplest example is a single scale cluster expan- 
sion but with so-called small/large field conditions [R2], in which one does not want to derive 
links between cubes of a large field connected component. Indeed large field regions are small 
for probabilistic reasons, but the vertices of perturbation theory (created by the action of 
functional derivatives such as those in (V.A.8)) may not be small there, hence it would be 
dangerous to blindly expand the clusters in the usual way. 

We want to study a theory with partition function such as 

Z(A) = Z ^ T ) > Z ^ T ) = [ d^ c {4>)xA4>)e- 9 ^ 1{4>)dx , (V.E.I) 
rcA •* 

where C is a propagator with decay at the unit scale, and 7(0) is some local interaction. 
Without describing a precise model we shall assume that the large field condition xt{4>) is 
such that the functional integral J (i / uc(0)xr(</ , )e ff J*A J ^ dx can be bounded by ATl A_r lc' r l 
where K is fixed and c is a small constant that can tend to zero with some coupling constant 
in 1 ((f))] one also assumes that the outcome of a functional derivative f b dx 5 ^ x ^ acting on 
gives a small factor but only in the "small field region" , hence if b £ V. 

The set I n is again the set of the cubes which pave our finite volume A. We fix a 
given subset V C I n and define, for / = {ij} G V n , ef = 1 if bi e T, bj G V and dist(b;, bj < M. 
M is some constant that will be fixed to a large value. Otherwise we put ef = 0. We define 
also rjf = 1 — efj . Hence ef = 1 means that bi and bj are "large field cubes" which are closer 
than distance M, and rjf = 1 means the contrary. 

The level two jungle formula will at the first level create connections whose clus- 
ters are automatically the "connected components" of V in the generalized sense that up 
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to distance M two cubes are considered connected. Then the second level will create ordi- 
nary connections with the propagator C. Therefore we introduce a first set of interpolation 
parameters {x}} and define 

H r ({x}})= U(xM + vf)- (V.E.2) 

iev n 

Remark that H T (1, 1) = 1, so that we can freely multiply Z(A, V) by H r (l, 1) without 
changing its value. We introduce the second set of parameters {xf} on the covariance C 
exactly as in (V.A.5). 

Then the formula (IV. 11) with m = 2 simply gives: 



z(a,t)= £ (n/^Ofn^fnw+^r 1 ^))) 



2— jungle 



J ^(h)(0) J] {/ dxd VXh (x) X b j {y)C{x, y) * * |xr(0)e 9 L^ dx 
l={ij}^S2—di 

(V.E.3) 

where d(ic T (h) is the normalized Gaussian measure with positive covariance 



C^(h)(x,y) = ^ tf) (h)C(x,y) (V.E.4) 

where y ^(h) is by definition 1 if x and y belong to the same cube, and is h t ' (h) if 
I = {ij}, x G bi and y G 6j or x G 6j and j/ G 6^. 

The only non-zero terms in this formula are those for which the clusters associated 
to the forest 3i are exactly the set of "connected components" r a of the large field region 
in the generalized sense (for M = it gives the ordinary connected components). Indeed 
they cannot be larger because of the factor rize;^ e /\ nor can they be smaller because of the 
factor riz^ {vf + e FfyT' 1 (h)) which is zero if there are some generalized neighbors (for which 
r]fj = 0) belonging to different clusters (for which hf-{h) = 0). Therefore the first forest in 
this formula simply automatically draws connecting trees of "neighbor links" connecting each 
such generalized connected component, but in a symmetric way without arbitrary choices. 
Remark that the factor ri/^j {vf + e \ hf "' 1 (h)) is bounded by one as expected for further 
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estimates. Then the second forest of the jungle automatically takes into account the first 
links, i.e. the existence of large field regions, to draw propagators connections, In this way 
the units connected by the full forest remain unit cubes instead of being either small field 
cubes or blocks of large field cubes*. For convergence of the thermodynamic limit one has 
then simply to check that all connections are summable (this is obvious for the finite range 
nature of the ef links), and that there is a small factor per link of the forest. For e links 
it comes from the probabilistic factor associated to the presence of the function xr and for 
ordinary links, it comes either from the functional derivative localized out of T, or from the 
large distance (at least M) crossed in the case of a link between two large field cubes. This 
distance induces a small factor through the decay of C. In conclusion there is no need to 
gain a small factor from functional derivatives localized in the large field region (which is 
usually impossible anyway), and the whole convergence becomes much more transparent, 
many combinatoric difficulties being hidden in the jungle formula itself. 

A concrete example in which this formula would somewhat simplify the argument is 
e.g. the single scale expansion of [L]; in [KMR] a three level jungle formula is used, in which 
the third forest hooks some cubes along the straight paths of the second forest propagators, 
in order to complete factorization in the context of a Peierls contour argument. 

F) Cluster or resolvent expansions with smooth localizations 

In some situations (for instance if momentum conservation is important) , it may be 
inconvenient to perform cluster expansions with sharp localization functions such as Xb m 
(V.A.8). But with smooth functions there is no naive factorization in (V.A.8). This is not 
a serious difficulty and it can be overcome e.g. by an auxiliary expansion (sometimes called 
a "painting expansion" ) on the interaction that creates protection belts around the clusters. 
But in this last section we remark that the Taylor forest formula (III. 2) also gives elegant 
solutions to this type of problems and treat the </> 4 interaction again as an example. 

Let 

1 = 5>6(*) (V.F.I) 
b 

* These blocks lead to unpleasant additional sums for where in the blocks functional 
derivatives really act, etc... 
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be a smooth partition of unity of IR by cubes of unit size. We assume that the support of Xb 
is contained in {x|dist(x, b) < 1/3}. The set of all 6's is then further restricted to the finite 
set Af of all cubes which are at distance zero of our finite volume A (hence include a unit 
corridor around it). Therefore 

Xa(x) = xa(x) J2 Xlix) (V.F.2) 
beAf 

so that the corresponding sums are finite from now on. 

Let us rewrite the (j) 4 theory of section V.A in terms of an intermediate ultralocal 

field: 



Z(A) = J dfici&e- 9 'L^ dx = J d/j CA (<P)du(V)e l ^L^^ v ^ dx (V.F.3) 

where dv is the normalized Gaussian measure with ultralocal covariance 5(x — y). 
Inserting the identity V.F.2 we get 

Z(A) = J d^ c (<P)dHV)e^L dx ^ b eJ^ 2 ^ v ^ (V.F.4) 

We define now the collection of Gaussian random variables {^(x)}, b G Af distributed ac- 
cording to the measure dfi({(f)b}) with covariance C(b,x ;b',y) = Xb(x)C(x,y)xb'(y), and 
the collection of Gaussian random variables {Vb(x)} distributed with degenerate covariance 
x ; 6', y) = 5(x — y). We have 

Z(A) = J ^({^})rfz/({H})e^/A dx ^^Ar^) 2 W^W (v.F.5) 

(to prove (V.F.5), remark that the right hand side of (V.F.4) and (V.F.5) are both Borel 
summable functions of g with identical perturbative series). We may now change slightly the 
boundary condition, dropping the restriction on the interaction range of integration. In other 
words Z(A) leads to the same thermodynamic variables than 

Z(A) = j ^({0 6 })^({H})e^^ i -^i' dxWi,)2(:B)n(a:) (V.F.6) 

Let us apply the Taylor forest formula (III. 2) to Z(A), interpolating both the co- 
variances C(b,x ;b',y) and F(b, x ;b',y), viewed as finite matrices with entries in Af whose 
elements are operators on L 2 (M d ). 
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It gives a result very similar to (V.A.7-8): 



2 ^= e n^( y o 



Yi,...,v„ i=i 



A ( y )= E (II f dh i) I ^c-(h)(0)^r-(h)(^) 

T on F ZeT ^° 

// 5 5 5 5 

dxdy I (a;)x6,- (y)C(^, y) 7-7T A , 7 ; + 5(x - y) — — 
V S(p bi (x) 5(p b .(y) 5V bi (x) 5 V bj (y) 

'v^J^beY I dx (<Pb) 2 (x)V b {x) p ^ 



l={ij}ET 



e 

where hi and bj are the two ends of the line /, and the sum is over trees T which connect 
together the set Y, hence have exactly \Y\ — 1 elements. The measures d\x C T (h)(0) an d 
du r T^(V) are normalized Gaussian measure on the restricted collections of Gaussian ran- 
dom variables {(p b (x)} and {V&(;r)} for b E Y. d^ C T^ has covariance Cy(h)(6, x ;b',y) = 
C{x,y)fiT{b, b'), where h T {b, b) = 1 and for b 7^ b\ h T {b, b') is the infimum of all parameters 
hi for / in the unique path in T joining the cube b to the cube b' . Similarly du T T ^(V) has 
covariance Ty (6, x ; y) = 5(x — y)hT(b, b'). 

Remark that the difference between (V.A.8) and (V.F.7) appears in the term 5(x — 
y) sv b S (x) sv b 5 .(y) wmcn would be zero if the functions Xb were sharp characteristic functions. 
This term corresponds to "emptying" the interaction on the borders of the clusters created 
so that full factorization occurs, even without sharp localization functions. 

Similar formulas can be written for other kind of interactions or for resolvent ex- 
pansions with smooth localization functions. We leave them to the reader. 
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